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",-J-". We propose a new method to study transverse flow effects in relativistic 
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q^ , nuclear collisions by Fourier analysis of the azimuthal distribution on an event 
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by-event basis in relatively narrow rapidity windows. The distributions of 
Fourier coefficients provide direct information on the magnitude and type 
of flow. Directivity and two dimensional sphericity tensor, widely used to 
analyze flow, emerge naturally in our approach, since they correspond to the 
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distributions of the first and second harmonic coefficients, respectively. The 
role of finite particle fluctuations and particle correlations is discussed. 
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In the study of specific processes such as flow phenomena, it is very important to find 
variables which are both easy to work with and have clear physical interpretations. Many 
different methods were proposed for the study of flow effects in relativistic nuclear collisions, 
of which the most commonly used are directivity [JT] and sphericity (three dimensional B 
or two dimensional ||) tensor. In the current paper we propose to use a more general 
method, which includes as a part both directivity and two dimensional sphericity methods 
in a natural way and gives a clear physical meaning to the analysis. Moreover, it provides 
ways to investigate complicated event shapes, which may not be described by an ellipsoid 
through three dimensional sphericity analysis. 

There exist two different approaches in flow analysis. The first one is to fit pt and/or 
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dN/dy distributions of different particles by assuming local thermal equilibrium and hydro- 
dynamical expansion. In this case some information about initial conditions, equation-of- 
state, expansion dynamics, freeze-out temperature can be extracted. (For recent analysis 
of such type of CERN SPS and BNL AGS data see 0,0). In spite of sometimes very 
good description of experimental data this approach has essential plague: it is very model 
dependent. Just the very existing of thermal equilibrium is still not well established. The 
alternative approach is to study the azimuthal event shapes (which are sensitive to flow) 
from experimental data without any model in mind. All the methods mentioned in the first 
paragraph including the method we are going to discuss in this paper are falling into this 
category. The information provided in this way can be used for comparison with different 
models: hydro-dynamical, statistical, cascade, multi-string model, etc. The study of a 
possible azimuthal anisotropy of the multiparticle production in high energy nucleus-nucleus 
collisions just has begun, and first indications of it have been seen at AGS 0, where the 
anisotropy of particle production in one region was noted to correlate with the reaction plane 
angle defined in another region. 

In our approach we investigate the shapes of azimuthal distributions, and try to detect 
and study of anisotropy of different kinds. The origin of the anisotropy could be different: 
hydro-dynamical flow due to pressure gradients, shadowing, both, etc. What they have in 
common is some collective behavior in the evolution of multiparticle production process. 
Consequently, we often use the term transverse collective flow for the phenomena we study, 
although we do not necessarily mean the hydro-dynamical flow. Our method is most useful 
for detecting an anisotropic collective flow. Though it can be used for a study of radial 
flow as well (see the short discussion in the Section |TJ) along with other methods such as an 
analysis of p t spectra. In order to detect complex three dimensional event shapes, we slice 
longitudinal variables, e.g., (pseudo)-rapidity, into different windows, and perform Fourier 
analysis in each window. Fourier coefficients of different harmonics reflect the different type 
of anisotropy (transverse collective flow). A three dimensional event shape can be obtained 
by correlating and combining the Fourier coefficients in different longitudinal windows. Some 
ideas of this approach have been employed to analyze Au on Au data of E877 collaboration 
at BNL AGS, when for the first time in this energy range transverse flow signals have been 



observed 0.|]. 

We present the Fourier expansion method in Section [TJ. We study the fluctuations due 
to finite particle multiplicity in section |T|. In section |V|, we address correlation between 
different windows and different harmonics. Finally, relations to other existing methods are 
discussed in section M. 

II. FOURIER EXPANSION OF AZIMUTHAL DISTRIBUTIONS 

The particle azimuthal distributions can be constructed from different quantities such 
as transverse momentum, multiplicity, or transverse energy in relatively narrow (pseudo)- 
rapidity windows. To get more detail of the event shapes one would like to divide the whole 
(pseudo)-rapidity acceptance into more windows. But the limitation is that there should be 
sufficient number of particles in each window. We will discuss this problem in next section. 
If needed, it is possible to introduce different weights for different particles, as is often done 
in the analysis of flow phenomena at lower energies. It is also possible to study specific 
functions of p t , multiplicity, etc; for example, if one wants to enhance the contribution from 
particles with high p t , one can study p1{4>) instead of Pt(<f>), or select particles with p t larger 
than certain value. In this way the method can also be used for a study of radial flow. To 
study flow effects in future RHIC experiments, where jet production can be important and 
will be a background to flow study, one can chose the range p t < c(pl e ), where constant 
c can be chosen from experiments. Or one can do just the opposite to study anisotropy 
due to jets. With particle identification, one can get the azimuthal function for certain 
type of particles, e.g., analyze protons and pions separately. In the following to be specific 
we generally discuss an azimuthal distribution of transverse momentum in certain rapidity 
window, but the procedure for other variables remains the same. 

We denote the azimuthal distribution function of the quantity under study with r(<p), i.e., 
dpx{4>) I 'd(j), where Pt{4>) is the total transverse momentum of particles emitted at azimuthal 
angle 0. The function r(0) can be constructed from experimental data event by event and 
be written in the form of Fourier expansion since r(0) is a periodical function: 
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The coefficients in the Fourier expansion of r(<p) are integrals of r(0) with weights propor- 
tional to cos(n0) or sin(n0) . For the case of a finite number of particles, the integrals become 
simple sums over particles found in the appropriate rapidity window: 
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x n = r((f))cos(n(j))d(j) = y2r u cos(n(f) u ), (2) 
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y n = r(0) sin(n0)d0 = Vr I/ sin(n0 I/ ), (3) 
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where v runs over all particles, and <p v is the azimuthal angle of the z/-th particle. Without 
any flow and neglecting fluctuations (we study the fluctuations in Section PTT[ ) the function 
r(4>) is constant, r(0) = x /(2tt) = l/(27r) Y^v r v All Fourier coefficients except a are zero. 
The appearance of transverse anisotropic flow in nuclear collision is expected from a non- 
zero value of an impact parameter. We define the longitudinal or beam direction as z-axis, 
and the transverse plane as x-y plane. The impact parameter vector in x-y plane (from the 
center of the target nucleus to the center of the beam nucleus) together with z-axis defines 
the reaction plane, and the reaction plane angle ip r (0 < ip r < 2ti) is the angle between 
x-axis and reaction plane. Each non-zero pair of the Fourier coefficients x n and y n give the 



value of non-zero component of the corresponding harmonic, defined as v n = \jx\ + y% : and 
the angle ip n (0 < ip n < 2ir/n) of the "n-th type" flow. 

x n = v n cos (ni> n ), (4) 

y n = v n sm(nip n ). (5) 

If flow exists, r(0) is no longer a constant, and the shape of the distribution is no longer 
a circle centered at zero (see Fig. [I]). Then the first harmonic coefficients correspond to an 
overall shift of the distribution in the transverse plane; such flow we call directed flow. The 
magnitude of flow is V\. With r(0) being the transverse momentum distribution, coefficients 
X\ and y\ are x and y components of the total transverse momentum of the particles produced 
in the window, while v\ is the magnitude of the total vector sum of transverse momenta. 
In the case of negligible fluctuations the direction of flow due to symmetry is to coincide 
with the reaction plane angle i\)\ = ip r (repulsive flow), or point in the opposite direction 
ipi = ip r + it (attractive flow). 



The non-zero second harmonic describes the eccentricity of an ellipse-like distribution. If 
one approximates the distribution by an ellipse as shown in Fig. [TJ, then the second coefficient 
v 2 carries information on the magnitude of the eccentricity. Quantitatively it is the difference 
between the major and the minor axis. The orientation of major axis ip2 (or, ip2 + fl", which 
give the same orientation for an ellipse) can be only ip r or ip r + tc/2. In the case i\)i = ip r , the 
major axis lies within the reaction plane; while ip2 — Vv + vr/2 indicates that the orientation 
of the ellipse is perpendicular to the reaction plane, which is the case for squeeze-out flow 
and may be expected in mid-rapidity window. 

Coefficient V3 could be non-zero for asymmetric nuclear collisions and represent the 
asymmetry of the flow due to different sizes of colliding nuclei, if the distribution can be 
approximately described as triangle-type distribution. In this case the symmetry dictates 
that the flow angle ^3 (or one of ^3 + 2mr/3, n — 1, 2) gets value ip r or ip r + ir/3. 

Forth - a rectangle-type deformation amplitude v 4 could be non-zero for rapidity windows 
close to the center of mass, where both squeeze-out, emitting preferentially more particles in 
the direction perpendicular to the reaction plane, and "side-splash" effects can be important. 
The flow angle ^ 4 would be if) r if the origin of the rectangle-type flow is the case as discussed 
above. But ip^ = ip r + 7r/4 can not be excluded from the naive symmetric reaction plane 
argument. 

We see that Fourier coefficients corresponding to different harmonics have a clear physical 
meaning. Perhaps the most interesting and unique feature of this approach would be to 
correlate flow effects between different types. But before trying to do that, we have to 
answer the question of how to deal with the finite multiplicity fluctuations, which have 
always been a big obstacle in flow analysis. Because it is difficult to resolve the reaction 
plane angle event by event due to fluctuations, all methods based on the determination of 
the reaction plane on an event-by-event basis remains rather ambiguous. It becomes very 
important to detect flow without using the information on the exact position of the reaction 
plane in each event. When consider the fluctuations, the value of v n becomes a distribution 
instead of a single number for events of the same characterization. Below we show that 
the distribution in v n is sensitive to anisotropic flow and can be used for flow detection. 
Information about the reaction plane position could be used if one considers the relative 



angular distance between the reaction plane and some other plane, also defined in the same 
event. For analysis of this kind we propose to study the correlations between reactions plane 
angles defined in different windows, or between reaction plane angles derived, possibly, from 
the same window but using different harmonics. 

III. FINITE MULTIPLICITY FLUCTUATIONS 

In this section we discuss problems arising in the analysis due to finite multiplicity 
fluctuations. Since two particle correlations are weak |J, below we treat all particles in the 
same event as being totally independent following certain distribution with variance ctq. If 
two particle correlations are non- negligible we expect larger fluctuations, but the general 
procedure remains valid. Although we study finite multiplicity fluctuations, we consider the 
case when the total number of particles in selected window is relatively large (N ^> 1). Under 
these assumptions, in the absence of flow the probability distribution of vectors (vi,if)i) is 
Gaussian from the central limit theorem: 

d 2 w 1 . v ? . 

ex P(-7rU (6) 
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with the peak of the distribution centered at the center of the coordinates. The variance of 
the distribution is: 

a 2 = al = al = N <r 2 >(cos 2 (^)> = ^- = iVa 2 , (7) 

where N is the number of particles in the window under study. 

For simplicity, we begin our study with flow of the first type, namely directed ("side- 
splash") flow. With the appearance of flow, each event can be characterized by the mag- 
nitude of flow and the flow angle. For events with the same parameters of flow the center 
of the distribution (||) will be shifted to the point characterized by (di,ip r ) in the polar 
coordinate system. It is reasonable to assume that in shifted distributions the statistical 
fluctuations (0) are the same. We have chosen here v\ as the parameter relevant to the 
magnitude of flow, which is directly related to the impact parameter of the collision. On 
the other hand, it is known from Monte Carlo studies that the impact parameter is strongly 
correlated with global variables used to measure the centrality of events (transverse energy, 
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multiplicity, zero degree energy, etc.). Using these correlations one can select events with 
the same magnitude of flow experimentally. Below we assume that the distribution in v\ is 
narrow and can be treated as a constant for events within the same centrality bins. 

The distribution in the flow angle ip% should be flat. That can be used experimentally 
to justify the equalization and correction applied to detectors. We can integrate over if)[ = 
■01 — VV; which gives: 

dw 1 f 27T ,n vf + v i ~ ZviVx cos(V4) 1 v\ + v\ viVi 
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where Iq is the modified Bessel function. The shape of the distribution (|8]) for a few values 
of parameters are shown in Fig. pj, where we plot dw/xdx using the notation x = V\/a and 
x = "Ui/cr. The extrema are defined by the solution of the equation 

V\ h(vivi/a 2 ) - vi I^vxvi/o 2 ) = 0. (9) 

It follows from Eq. ^| that the distribution (Eq. [8D has a local minimum at V\ — if V\jo > 
v2. Thus the observation of such a minimum in the data is a sufficient condition for the 
appearance of flow. But such a minimum is not required for the extraction of V\ and a by 
fitting Eq. flD to data. 

It will be shown in Section |VC| that in the case of the analysis of azimuthal pt distri- 
butions v\ equals N{p x ), where (p x ) is the mean transverse momentum gained by a particle 
due to transverse flow (x-z plane is the reaction plane here). In this case 

vi N(p x ) ( Px )V2N 



(10) 

We see from Eq.([T0|) that, one has to slice the rapidity windows experimentally to accept 
reasonably large multiplicity (JV) in each window to make sure that v\ja is not too small, 
and an available statistics allows to perform a reliable fit. The "relatively narrow" rapidity 
windows discussed in the text are precisely chosen this way. 

Applying the same consideration to the distribution of the coefficients corresponding to 
higher harmonics, one can see that they should have the same behavior (Eqs. 0— 1|) but with 
their own parameters Hi. Note, that because of (cos 2 (0)) = (cos 2 (n0)) one has the same 
values of a for all harmonics. 



The functional form (|8D, which allows to have a local minimum of the distribution at zero, 
was proposed almost a decade ago for the distribution of flow polar angles extracted from 
the three dimensional sphericity matrix [HJ. For the two dimensional transverse sphericity 



matrix a comprehensive study was done in Ref. P,pJ|. Here we want to emphasize that the 
same functional behavior with the same parameter responsible for finite multiplicity fluctu- 
ations is valid for the distribution of Fourier coefficients of any harmonic order. Therefore 
one would expect more stable and sensible results fitting simultaneously all dw/v n dv n dis- 
tributions with free parameters v n (n = 1, 2, 3, 4, ...) and the same parameter a. 

In order to be able to compare the magnitudes of flow in different rapidity windows or 
centrality bins, or to compare different experiments, it is useful to consider the values of 
v n normalized to v = (xq) = N(p t ), that is v n /vQ. The ratio V\/vq has a clear physical 
meaning as the ratio of transverse momentum gained by a particle due to transverse flow to 
the mean transverse momentum (x-z plane being the reaction plane): 

?i = H (id 

IV. CORRELATIONS 

The parameter v n and a can be extracted from experimental data by fitting the distri- 
butions dw/v n dv n by the functional form of Eq. || The non-zero v n reflects the existence 
and magnitude of flow of "n-th type" . It is interesting to study whether the detected flow 
in one rapidity window is correlated with flow in a different rapidity window, or whether 
the flow of one kind of particles is correlated with the flow of another particle type. The 
study of some of such kind of correlations was proposed as proof of an existence of flow 
in Ref. [1] (where the authors discuss in our notation the correlations of ip\ from different 
rapidity windows) and in Ref. (T^j (correlations of ipz), in our more general approach along 
with these correlations we can study very interesting correlations between flow of different 
types defining the overall flow picture. In general, we can study the correlations between 
any rapidity windows and any harmonics, (v' m ,mip' m ) and (v n ,nip n ). Denoting the relative 
angle by <fi = ip' — ip n and averaging over the reaction plane angle we have the distribution: 
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The observation of correlation expressed in Eq. 12 could be very helpful for understand- 
ing the nature of flow and get a global picture of three dimensional event shapes. The 
correlations between flow of different types in the same rapidity window can be biased by 
auto-correlations due to finite multiplicity To avoid the problem one can study the cor- 
relations using as a reference a flow angle defined in another rapidity window. The other 
way to suppress these pseudo-correlations could be to calculate higher harmonic term after 
removal of lower harmonic terms event-by-event, e.g., to perform a shift of the center of 
particle distribution to the point (vi,ipi) and then calculate (t> 2 , 2t/> 2 )- 

The correlation between different rapidity windows can be numerically calculated using 
the fit parameters to the distributions (|8]). In this case a comparison with experimental data 
will be a check of the consistency of the method. Note that in spite of rather complicated 
form (0), the real calculation can be easily done numerically. Some analytical calcula- 
tions for specific parameters can be done using the technique presented in the Appendix to 
Ref. flTTl. 



V. RELATION TO OTHER METHODS 

A. Directivity and sphericity 

The vector Q, which is often referred as directivity, was originally introduced in the flow 
analysis by Danielewicz and Odyniec Q: 

Q = J2 UJ »Vt;v, (13) 

v=l 

where v is particle index and uj v is a weight. It coincides with the vector vx = (vi, ipi) in the 
Fourier analysis by taking io v — 1. One of the results of Ref. |1[ concerning the width of the 
distribution in Q is directly related to the modification of the original Gaussian distribution 
due to flow, which we discussed above. 



The relation of the second harmonic coefficients to the usual parameterization of two 
dimensional sphericity tensor is also straightforward. Two dimensional transverse sphericity 
tensor, which is relevant to our study, has three independent parameters, and usually is 
written in the form: 

^ S 1 1 + e cos 29 e sin 29 \ 

Sii = 2^u v p i . v p j . tV = -\ , (14) 

v=\ z y esin20 l-ecos2#y 

where i, j = x,y. If we chose u v = r u /(p t - u ) 2 the sphericity tensor becomes 

N 

S ij = Yl T " COS 0",« COS ^,h ( 15 ) 

v=\ 

where <p v ^ is an angle between the momentum of v-th. particle and i-th axis. In this case 
S = trSij = J2v r v = vq. By computing the eigenvalues one can derive that the eccentricity 
e = v 2/v and thus is directly connected to the magnitude of the second harmonic coefficient 
in Fourier analysis. 

B. Reaction plane determination 

Experimentally, to study the flow in the reaction plane (so called side-splash flow), tra- 
ditionally a directivity (an angle ip{) is calculated from the whole available acceptance to 
define the reaction plane. To define the accuracy of the procedure one should consider the 
distribution in an angle ip which is defined as ip — ip r — tpi : 

dw f vdv , v\ + v\ — 2v\V\ cos ib . 

- " exp(- '' 



dip J 2na 2 ^ v 2a 2 ' 

1 v\ wicos^v/tt i) 2 sin 2 ^ v x cos$ 

= 27 [eXp( -2^ } + aV2 eXp{ —2^ ){1 + erf( T7T ))] - (16) 

As long as the parameters of flow are determined by fitting the experimental distribution 
to the form ([8]), Eq.(|T6|) can be used to answer the question of how close the angle ipi is to 
the true reaction plane angle ip r . It can be generalized to n-th type flow angle dispersion by 
replacing v\ with v n , and ipi with mp n . 
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C. Mean p x 

Any transverse flow would influence the pt distribution of observed particles. Theoretical 
calculations |12| usually present a mean projection of the particle transverse momentum into 
the reaction plane (p x ) as a function of rapidity (it is assumed here that x-z plane is the 
reaction plane). The projection can be non-zero due to transverse flow. Generally, as 
discussed above, the directivity is used to define the reaction plane, and then in different 
rapidity bins transverse momenta are projected into this plane. (One should be careful here 
to avoid autocorrelation. For details see Ref. JI|.) Below we use the notation p x > for such a 
projection instead of p x to note that it is not the projection into true reaction plane. Due to 
the fluctuations in ijj\ around the true reaction plane angle (p x >) does not equal (p x ). In our 
approach, we get our fitting parameter V\ without determining of the reaction plane on an 
event-by-event basis. Taking into account that vi = ^ Pt;u, one can see that v\ is exactly 
the true (p x ), multiplied by the total number of particles in the rapidity region. 



lr ^ (vi - V 2 

o — 2 / dvii>icos(vi,vi)exp( r-j 

liter* J la z 



(p x ) N = ^—^ I d\\V\ cos(vi, vi) exp( — ) = v 1 . (17) 



The relation between (p x ) and (p x i) is very simple [jl| (p x ) = {p x ') {cos(ijj xx i)) , where ipxx' 
is the difference between true and defined reaction plane angles ip xx i = ip r — ipi. In our 
approach 



1 f — ~ 

(cos(^O) = 2 — 2 / dviCOs(vi,Vi)exp( 



2a 2 



^i -- t ^i \ d ,o /o. o. vf 



eM~)iFi(3/2;2;-±), (18) 

y/2n <J 2a z 2o A 

where ii 7 \(a;/3;7) is a confluent hypergeometric function. The behavior of (cos(Vw)) as a 
function of v\ja is shown in Fig.|^. 

D. Azimuthal two-particle correlations 

An azimuthal two-particle correlation function C(ip) = Peoria) / Puncori^) , where ip is 



the angle between two particles, was proposed in Refs. |113l , |14| to study transverse flow 
as a method which is free from uncertainty in the reaction plane determination. Here 
P C or(4>) is observed pair distribution and P U ncor(4>) is the distribution generated by "event 
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mixing" . The formalism used above for Fourier analysis permits us to estimate the azimuthal 
two-particle correlations in the same way as it was done for the distributions of Fourier 
coefficients. The origin of azimuthal correlations lies in the overall shift of the particle 
distribution in p t -space due to flow. The parameter which controls the magnitude of the 



correlations is (j> x ) / y (j>t) ■ The correlations increase with increasing values of this parameter. 
Note that the corresponding parameter (which defines the magnitude of the effect) in our 



approach is V\ja = N{p x )(y)/JN{pf)/2, which is \/N times larger. The difference, as 
always, originates from the fact that in the study of collective phenomena the fluctuations 
in the collective variables are less than the fluctuations in variables describing individual 
particles. 

If one assumes a Gaussian form for the one-particle pt distribution, the expression for 
C(ip) coincides with the expression for the correlation between two flow angles defined in 
different regions, after a substitution of corresponding parameters. The statement made 
above about finite multiplicity fluctuations becomes explicit in this case. 

VI. CONCLUSION 

To summarize we propose to study transverse anisotropic collective flow by analysis of 
azimuthal distributions in different rapidity regions. We show that Fourier expansion of the 
azimuthal distributions is an appropriate tool for such an analysis. It is shown that the 



distributions in v n = \jx\ + y^, where x n and y n are n— th harmonic Fourier coefficients, 
are sensitive to transverse flow. The classification of flow of different types leading to a 
deformation of azimuthal distribution corresponding to different Fourier harmonics appears 
naturally in this method. The important feature of the approach is that it is free from uncer- 
tainties in a determination of the reaction plane on an event-by-event basis. The dispersion 
of the flow angles is derived through the discussion of the finite multiplicity fluctuations. 
Correlations in the flow angles determined for flow of different types are proposed as a tech- 
nique for flow analysis. It is shown that many existing methods are within our framework, 
either the first or the second harmonics. But our approach provide additional information, 
from which a three dimensional event shape can be established. 
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FIGURE CAPTIONS 

1. Function r(4>) (solid line). Approximation by ellipse (dashed line). 

2. The distribution dw/xdx for different values of x. The curves marked by 1,2,3, and 4 
correspond to x — 0.5, 1.0, 1.5, 2.0, respectively. 

3. (cos((f> xx i)) as a function oiv\/a s 
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